Can entropy save bacteria? 
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This article presents a physical biology approach to understanding organization and segregation 
of bacterial chromosomes. The author uses a "piston" analogy for bacterial chromosomes in 
a cell, which leads to a phase diagram for the organization of two athermal chains confined 
in a closed geometry characterized by two length scales (length and width). When applied to 
rod-shaped bacteria such as Escherichia coli, this phase diagram predicts that, despite strong 
confinement, duplicated chromosomes will demix, i.e., there exists a primordial physical driving 
force for chromosome segregation. The author discusses segregation of duplicating chromosomes 
using the concentric-shell model, which predicts that newly synthesized DNA will be found in 
the periphery of the chromosome during replication. In contrast to chromosomes, these results 
suggest that most plasmids will be randomly distributed inside the cell because of their small sizes. 
An active partitioning system is therefore required for accurate segregation of low-copy number 
plasmids. Implications of these results are also sketched, e.g., on the role of proteins, segregation 
mechanisms for bacteria of diverse shapes, cell cycle of an artificial cell, and evolution. 
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I. INTRODUCTION 

Chromosomes are a cornerstone of fundamental pro- 
cesses of any cell, and harmony between their physi- 
cal properties and biological functions is an evolution- 
ary consequence. In bacteria, two defining processes of 
the cell cycle, DNA replication and chromosome segre- 
gation, progress hand-in-hand. The model system Es- 
cherichia coli contains a single circular chromosome in a 
rod-shaped cell. DNA replication starts at a unique ori- 
gin of replication (ori), creating a replication bubble that 
grows bidirectionally, and the two replication forks meet 
at the terminus (ter) located at approximately the op- 
posite clock position of ori on the circular chromosome 
(Fig. [TJ. In slowly growing cells, there is a one-to-one 
correspondence between each complete round of replica- 
tion and the cell cycle. In fast growing cells, on the other 
hand, the cell divides more frequently than the progres- 
sion of the forks from ori to ter and, thus, new replica- 
tion has to initiate before the completion of the previous 
round of duplication, leading to "multifork" replication. 
In either case, replication and segregation are coordi- 
nated with the growth of the cell, and recent visualization 
experiments have begun to reveal how replicating chro- 
mosomes move and segregate during the bacterial cell 
cycle. Although details vary from organism to organism, 
the observation common to all bacteria of rod-shaped 
cells studied so far shows directed movement of duplicat- 
ing chrom osomes along the long axis of the cell [as seen 
in E. coli ([Bates and Klecknerj |2005l |Elmore et a?.[|2005' 



Nielsen et al. 2006 Wang et al. 


|2006|), Caulobacter cres- 


centus (Viollier et al. 2004), Bacillus subtilis (Berkmen 


and Grossman 


2006) and Vibrio cholerae (Fiebig et al. 


2006 Fogel and Waldor| 2006)] 


Expcrimental data also 
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suggest that, once duplicated, there is a fair degree of 
correlation between the relative clock position of circular 
chromosome and their average positions along the long 
axis of the cell, which often (but not always) shows the 
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FIG. 1 Schematic illustration of the bacterial cell cycle. 
Replication begins at the 12' position (ori) on the circular 
chromosome. Replication forks grow bidirectionally and meet 
at the opposite clock position (ter). During the bacterial cell 
cycle, replication and segregation progress hand-in-hand, cou- 
pled with the growth of the cell. 



principal linear ordering of chromosome (Teleman et al. 



1998 |Viollier et al. 2004 Wang et al. 20061 



Perhaps the most influential model on chromosome 
segregation s o far was proposed by |Jacob, Brenner, and| 
Cuzin ( 1963 1 in their seminal paper on the replicon model 
of E. coli: if replicating chromosomes (e.g., duplicated 
ori's) are attached to the elongating cell- wall membrane, 
they can be segregated passively by insertion of mem- 
brane material between the attachment points during 
growth. The model is intuitive and elegant - but wrong. 
(See Sec.jlLFj). 

Since the work of Jacob and colleagues, two classes 
of models have been proposed to explain the observed 
directed drift of duplicating chromosomes: (a) biological 



mechanisms such as DNA replication (Lemon and Gross- 
man] 20011, co-transcriptional translation of membrane 



proteins ( Woldringh[ |2002[ ) , RN A transcription ( jDworkin 
and Losick 2002 [ ); and (b) physical/mechanical driv- 



ing forces such as mechanical pushing between chro- 



mosomes (Bates and Kleckner 



2005 



Kleckner et al. 



2004 ) or conformational entropy of duplicating chromo- 



somes (Danchin et al. 2000 Fan et al. 2007 Jun and 



Muldei 



20061. These two classes of models are not mu- 



tually exclusive. 

The second class of model (b) is important in that it 
provides a basic physical framework to examine the ne- 
cessity and the roles of proteins or machinery (if any) 
involved in the cellular process and, if so, what condi- 
tions these putative biological players need to satisfy in 
order to successfully carry out their biological functions. 
For instance, when the cell is lysed, the bacterial chromo- 
some expands to several times the size of the encapsulat- 
ing cell ( |Cunha et aT||2001| (Fig.|3}\). This observation 
is puzzling since duplicated chromosomes must stay seg- 
regated without mixing until cell division for the viability 
of the cell, and yet the natural size of each chromosome 



is larger than the confining volume of the cell itself; why 
do the chromosomes not mix when they should occupy 
as much volume as possible? Thus, if the basic physi- 
cal principles predict that the duplicated chromosomes 
should mix inside bacterial cell, this would explain why 
putative proteins must be actively involved in segrega- 
tion of replicating chromosomes. However, if the physi- 
cal analysis predicts otherwise, this would imply the ex- 
istence of a primordial physical driving force underlying 
chromosome segregation in bacteria and that active bi- 
ological mechanisms may be present for other reasons 
(e.g., when the cell cycle is perturbed by external stress). 

The purpose of this paper is to present a compre- 
hensive and rigorous physical model of bacterial chro- 
mosomes that makes experimentally testable predictions 
based on measurable physical parameters of the bacte- 
rial cell and its chromosomes. 1 We have developed our 
model with the idea that physical properties of the chro- 
mosomes must interplay with the shape and size of the 
cell, and bacteria might have evolved to learn how to 
regulate these parameters to ensure proper segregation 
of chromosomes. 

This article is organized as follows: First, we begin 
with the concept of entropy, providing specific examples 
to explain how order can emerge from disorder driven 
by entropy. Second, we explain, using a piston analogy 
of bacterial chromosomes confined in a cell, how the 
geometric conditions of the cell interplay with the 
physical state of the confined chains. Based on this 
analogy, we present a phase diagram for the organization 
of self- avoiding chains in a box whose dimensions are 
defined by two lengthscales (width and length). To 
our knowledge, this is the first direct application of 
the de Genn es-Pincus blob theory ( de Gennes) 1979 



Pincus 1976) to closed geometries with more than one 
length scale (as opposed to, e.g., open channels, slits, or 
spheres, which can be characterized by a single length 
scale). 2 Then, we critically examine the model system 
E. coli and predict that duplicated chromosomes of 
this organism should segregate spontaneously (with the 
ordering of individual chromosomes being principally 
linear unless interfered by tethering to the wall). We 
also predict the opposite behavior for plasmids (namely, 
random distribution) because of their small sizes. 
We discuss various implications concerning entropy- 
driven organization and segregation of chromosomes 
in bacteria of diverse shapes as well as in an artificial cell. 



1 For physical approach to eukaryotic chromosomes, the reader is 
encouraged read iMarko and Siggial (119971 . 

2 In his celebrated book, Healing Concepts in Polymer Physics, 
de Gennes noted "This analysis can be extended to chains that 
are squeezed in slits and to other geometries provided that the 
confi ning object is cha racterized by a single length D." [see p. 
51 inlde Gennesl rtl9T9b] . 
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FIG. 2 Subtle aspect of entropy. (Molecular dynamics sim- 
ulation.) (a) Two species of particles (iVi = N2 = 200), 
initially separated by a wall, mix as the wall is removed, (b) 
If chain connectivity is introduced to this system by linearly 
connecting the particles of the same species, the two chains 
segregate. (iVi = N 2 = 512) 



II. PHYSICAL BIOLOGY APPROACH TO BACTERIAL 
CHROMOSOMES IN A CELL. 



There are also other examples in soft condensed mat- 
ter physics where entropy leads to ordering. For instance, 
Onsager's hard-rod model (19491 of the nematic-isotropic 
transition of liquid crystals is based on a similar physical 
insight of the trade-off between loss of orientational en- 
tropy and gain in positional entropy when hard rods are 
oriented parallel to one another. Crystallization of hard 
spheres (Alder and Wainwright 1957 Wood and Jacob- 



son 



1957) is another example where ordering allows a 
larger room for fluctuations. 

The above examples clarify what entropy really mea- 
sures, namely, the degrees of freedom, or the size of the 
phase space, of the system. For the main theme of this ar- 
ticle, the specific case in Fig.[2]points out the importance 
of chain connectivity, a keyword in polymer physics. In 
this example [Fig.[2jb)], chain connectivity defines a con- 
formational space f^B within the configuration space £Ia, 
and within £!b the typical conformations of the chains 
are those de- mixed. This emergence of order from disor- 
der due to entropy is our starting point of understanding 
chromosome segregation in bacteria. 3 



A. Entropy measures the degrees of freedom of the system 

In his influential book, What is life?, |Schrodinger| 



( 1944 ) asserted that the only thcrmodynamically equi- 
librium state of a living being is death. Defining entropy 
as a measure of 'disorder', in the last chapter 'Is Life 
Based on the Laws of Physics?', he struggled to explain 
how order (life) can be achieved from disorder. This as- 
sociation of life with minimal entropy, however, can be 
misleading. 

In Fig. |2^a), we show snapshots from molecular dy- 
namics simulations, where two species of Ni and N 2 par- 
ticles in equilibrium (blue and red, respectively) in a long 
rectangular box, initially separated by a wall. As we re- 
move the wall from the system, the two species start to 
mix. The driving force of this process is the well-known 
"entropy of mixing," which can be estimated as 



AS A = fc B (lnfi A -lnfi A> o) 



iViln 



N 



AT? In 



N2 
N 



(1) 

> 0, (2) 



where ks is the Boltzman constant, N = Ni + N2, and 
W A = N\/N 1 \N 2 \ (Qa,o = 1) denotes the total number of 
configurations of the system after (before) mixing. 

Entropy, however, is more subtle than a simple mea- 
sure of disorder. To see this, let us start with the mixed 
state of Fig. |2ja) and connect the particles of the same 
species, creating two long linear chains, one painted with 
blue and the other with red [Fig. |2^b), left]. Another 
important condition is the excluded-volume interaction 
between the particles. The reader is encouraged to per- 
form this simple computer simulation, and she or he will 
see that the two chains de-mix, i.e., that "order" emerges 
out of disorder (for the exact condition of segregation, see 



B. Piston analogy of bacterial cells and chromosomes. 

Imagine a cylinder containing two long linear molecules 
with excluded-volume, closed by two pistons which are 
initially far apart and do not perturb the chains, as shown 
in Fig. |3jh) top. In this dilute, confined solution of poly- 
mers, the chain conformations are well-described by the 
blob model; each chain is a series of blobs, stretched along 
the long axis of the cylinder without large-scale backfold- 

The origin of this 



de Gennes 


1979 


Pincus 


1976 



lapping blobs costs free energy ( Grosberg et al. 1982 Jun 



et al. 2007 1 . As we shall discuss below, it can explain 
the principal linear chromosome organization observed in 
rod-shaped bacteria. 

Here, the notion of "blob" is the central concept of 
our approach to bacterial chromosomes. It can formally 
be defined as an imaginary sphere within which the local 
chain segment densities are correlated ( |Doi and Edwards 
19861. In the above-mentioned case of dilute solution in 



a cylinder, the blob size £ is defined by the width of the 
cylinder D. It is important to note that a blob is also a 
unit of free energy of the chain and each blob contributes 
~ 1 fcffT to the free energy, regardless of its size (Jun 
et al. . 20071. As we shall discuss in Sec. Ill for bacterial 



chromosomes, the closest concept to the polymer blobs 
in the literature is macrodomains (Espeli and Boccard 



2006[ ) or "topological domains" ( |Deng et aL[|2005[ ). For 



3 Note that polymer physics was still in its infant stage during 
Schrodinger's time. For instance, Flory's magnum opus, Prin- 
ciples of Polymer Chemistry jl953| ), appeared almost 10 years 
after What is Life? Il944l. 
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(a) 



(b) linear ordering of chain conformation 




5 "blob" 



pushing -^--^ 



mixing 



random ordering of chain conformation 

FIG. 3 Piston analogy of bacterial chromosomes in confine- 
ment, (a) When the cell is lysed, the isolated chromosomes 
expands to several times the size of the cell. For our analy- 
sis, we used the published values of the nucleoid size b efore 

" [19991 



and after lysis, 0.24 /im (Woldringh and Odijk 



3.3/im (Romantsov et al. 
strain 



and 



2007 1, respectively, for E. coli B/r 
(b) Two athermal chains (i.e., with excluded-volume 
in good solvent) are confined in a cylinder whose pore width 
much smaller than the size of an unperturbed chain. The 
organization of the chain depends on the aspect ratio of the 
cylinder and the density of the monomers. 



instance, the latter denotes the largest independent chro- 
mosomal subunits, whose topological changes (such as 
cutting) do not affect other neighboring domains of the 
chromosome (for example, one may imagine loops that 
are independently anchored to one another). Indeed, it 
is also our view that chromosomes are best understood as 
a series of blobs, where each blob acts as an independent 



structural unit [see, also, Romantsov et al. (2007)] 



If we compress the pistons further, the two chains will 
touch one another and then partly overlap. This, how- 
ever, is not a favorable situation because overlapping 



chains repel each other 4 ( 


Bates and Kleckner 


2005 Jun 


and Mulder 


2006 


|Kleckner et al. 


2004 


) anc 


retract to 



occupy each half of the cylinder. Since the total acces- 
sible volume within the cylinder has decreased by the 
moving pistons, the total internal energy has to increase. 
In other words, the blob size will decreases accordingly 
(see above) [Fig. ^b) upper box] . 

We emphasize that the principal organization of 
the retracted chains is still linear, and this average 
linear ordering is maintained as long as the chains stay 



segregated. If the piston continues to move and increase 
the compression of the chains, the blobs also continue 
to become smaller to reflect the increased free energy. 
Once the blobs reach their critical size, the chains gain 
more conformational entropy by mixing with each other, 
and the chains lose their principal linear ordering and 
instead become a random walk of the connected blobs 
(Fig. [3jb) lower box). 



C. Phase diagram for mixing and de-mixing of polymer 
solution in a closed geometry. 

The piston analogy illustrates how the physical prop- 
erties of the chromosomes interplay with the confining 
geometry and its volume. We have obtained by rigorous 
analysis the phase diagram in Fig. [4] which describes the 
organizations of two chains in a rectangular box charac- 
terized by two length scales, width D and length L of the 
box. 

The dimcnsionless quantity x — Rf /D is the ratio be- 
tween the size of the unperturbed chain (Rf) and the 
width of the cylinder; thus, the x-axis represents how 
strong the confinement is. The y-axis indicates how dense 
the polymer solution is in the box, because the correla- 
tion length (blob size), £tmik, decreases monotonically as 
a function of the monomer density of the polymer in so- 
lution (see Appendix [A| . 

Our phase diagram is divided into two regions, one 
where chains segregate (blue), the other where they mix 
(red) (Fig. [4]). The diagram is further divided into six 
regimes. The boundary between segregation and mixing 
is given by the following geometric condition 



k 



L 
D 



D 

£bull 



(3) 



i.e., this is where the aspect ratio of the box, k, equals 
the ratio between the width of the box and the size of 
the blobs. The rest of the phase diagram can be best 
understood as we move along the horizontal dotted line 
in Fig. [4] Since £buik depends only on the chain density, 
for chains with fixed contour length, this line represents 
the change of the aspect ratio of the confining box (k) 
from filamentous (k > 1) to slab-like (k < 1) through 
cubic/spherical (k — 1) shapes, keeping constant the vol- 
ume of the box. Note that the longer the box is the better 
the chains segregate, and that the chains can readily mix 
m a sphere (k = 1) flJun et al.||2007| ). 

We leave a more complete description and technical 
details of the phase diagram to Appendix [A] 



4 This chain demixing can be understood intuitively. If we freeze 
the red chain, its inner envelope volume can be considered as a 
system of fixed obstacles. Thus, the blue chain moves away from 
the red chain, because the blue chain's gain in conformational 
entropy exceeds its entropic loss by occupying smaller volume 
free of the obstacles. This phenomenon is due to the chain con- 
nectivity. 



D. Application to bacterial chromosomes. 

Using the phase diagram presented in Fig. [4j we 
can predict whether a specific organism can segregate 
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FIG. 4 Phase diagram of a 2-chain polymer solution confined in a closed geometry of dimensions DxDxL. The x-axis represents 
the shape of the geometry. The y-axis represents the concentration of the polymer solution. Regimes with a red background 
represent mixed states, while those with a blue background represent segregated states. The dotted horizontal line shows 
the effect of aspect ratio of confinement on two chains of constant length N. Shown along this line are a few representative 
chain conformations (b)-(g). For cubic confinement with side D = Rf (a), two chains mix at the free-energy cost of order 
k B T ( jGrosberg et aZ. | [1982) |Jun et"aT) |2"007) . (h) Same as (d), but the two chains are in partial overlap of penetration depth a. 
Due to compression, the correlation length £buik of the confined chains is smaller than the width of the rectangular box D (see 
the stacking blobs). [Ro is the equilibrium end-to-end distance of an unperturbed chain in a long cylinder, and k the aspect 
ratio of the box.] See Appendix |A| for more technical details. 



its chromosomes using only conformational entropy - 
with no active control mechanism. To this end, we 
examined the model system E. coli B/r (H266), one 
of the most well-studied and well-documented bacterial 



strains (Woldringh and Odijk 1999). When in steady 
state, slowly growing (doubling time of 150 minutes) E. 
coli of this strain have a rod-shaped cell with hemispher- 
ical caps, with average length is 2.5/mn and the width 
0.5/im. However, the actual dimension in which chro- 
mosomes are confined inside a cell, namely the nucleoid, 
is smaller because of the well-known phenomena of nu- 
cleoid compaction (Woldringh and Odijk 1999| ). As a 
result, these values are estimated to be D = 0.24/um and 
La = 1.39/im for new-born cells containing a single chro- 
mosome 5 . 

On the other hand, the correlation length £buik is a 
more difficult quantity to measure experimentally, be- 
cause a number of factors contribute to chromosome or- 



ganization in bacteria (e.g., molecular crowding, super- 
coiling and DNA-protein interactions) (Stavans and Op- 



penheim 2006). Nevertheless, we can safely assume that 
the lower bound of £buik should be the persistence length 
of dsDNA, t p = 50nm. 6 Note that recent experimen- 
tal studies have shown that there are as many as ~ 400 
topological domains in E. coli chromosome, correspond- 
ing ~ lOkb per domain ( Deng et al. 2005 1 . If we assume 
that each domain occupies d 6 of a nucleoid volume, we 
obtain d ~ 53nm, about the same as £ p of bare DNA. 

Very recently, Krichevsky and colleagues have mea- 
sured £buik of nucleoid isolated from the same E. coli 
strain [B/r (H266)] using a more direct method of 



fluorescence correlation spectroscopy (FCS) (Romantsov 



et al. 2007[ ). Their careful analysis of the amplitude of 



the FCS correlation functions of randomly fluorescent 
labeled nucleoids revealed «50kb of structural unit 
which corresponds to d ~ 87nm of the B/r str ains in vivo 
[based on the nucleoid volume «0.06/^m 3 (Woldringh 



' The population average is L = 1.8um ([Woldringh and Odijk[ 
|1999| with L « 0.7L jKubitschek] 11981) , and we assumed a 
linear relationship between cell size and nucleoid length. 



For supercoiled DNA, the persistence l ength of the plectoneme i s 
estimated to be twice that of dsDNA (Marko and Siggia 1995 i. 



6 



R f 1000- 



Sbulk 



100- 



10- 



VI 



plasmid 




(chromo some) 



0. 



10 



100 



1000 



D 



FIG. 5 Segregation of DNA in E. coli. E. coli is in the 
compressed-cigar regime III (the bars represent 50% range 
of the average values x — 13.6 and y = 37.5, which we added 
to show the robustness of our estimate). Thus, the duplicated 
chromosomes have principal linear organization in the cell and 
segregate spontaneously to maximize their conformational en- 
tropy. On the other hand, plasmids are typically orders of 
magnitude smaller than the chromosomes, i.e., the point for 
chromosome should move on the phase diagram parallel to 
the y — x line and enters the mixing regime (e.g., V). Thus, 
during the course of evolution, plasmids may have acquired an 
additional, more active segregation mechanism other than en- 
tropic repulsion to ensure each daughter cell receives at least 
one copy. 



small "chromosomes" from a physical point of view (note: 
typical chromosomes are ^10 3 kb). At the insertion of 
two such plasmids in bacteria, the total amount of DNA 
(and, thus, £buik) remains practically constant, although 
the total number of chains increases to n = 4. In our 
phase diagram, this corresponds to changing both x and 
y of chromosome by a factor of ~ 10 1 (since Rp ~ N 3 ^ 5 ), 
i.e., the plasmids enter the mixing regime V 7 and they 
will distribute randomly inside the cell. Indeed, recent 
experimental results on the mobility and distribution of 
synthetic minimalizcd RK2 plasmid lacking the partition- 
ing system are fully consistent with our prediction of ran- 
dom distributions (iDerman et al. 2008). 



The quantitative reasoning above has broader im- 
plications for segregation strategies and copy-number 
control of plamids ( |Nordstrom and Gerdes 2003}: (i) An 
active segregation mechanism, especially, for low-copy 



number plasmids 1 


Adachi et al. 


2008 


Gerdes et al. 


2004 


)• or (ii) 



2006 Dcrman et al. 



i.e., without a protein-based segregation mechanism, 
plasmids may produce multiple copies so that, despite 
their small sizes, their chance of survival increases 
upon cell division. In random segregation, copy-number 
control is an interesting and important issue [See, for 



example, Brenner and Tomizawa| ( 1991 ) and references 
therein! 



and Odijk 1999)]. They also calculated the physical 
size of the unit based on less direct measurement of the 
diffusion constant of isolated nucleoids. Their estimated 
value was d — 70 ± 20nm, again, in good agreement 
with d = 87nm above. Indeed, we have used these 
parameters including the average size of the isolated 
nucleoids Rp = 3.3/xm (Romantsov et al. 2007[ ) (from 
the average spherical volume of natively supercoiled 
nucleoid, 18 /urn 3 ) and obtained x = 13.6 and y — 37.5 
- the E. coli B/r strains are in the segregation regime 
II (Fig. [5}. Thus, the replicating chromosomes could 
segregate purely driven by conformational entropy, and 



they will remain demixed (Bates and Kleckner 2005 



Jun and Mulder 2006) 



F. Segregation of replicating nucleoid. 

So far, we have discussed the physical conditions in 
which duplicated chromosomes will stay segregate, i.e., 
the conditions for a primordial driving force for chromo- 
some segregation. How then can we explain observed 
trajectories? 

Figures. [6] and [7] show experimental data of the aver- 
age positions of the chromosome loci in E. coli (Bates and 
Kleckner] [2005] and C. crescentus ( | Viollier et al.\ |2004[ ) . 
In E. coli, for instance, ori first moves towards the mid- 
cell position and, then, replication starts. The duplicated 
ori's split and move, 8 on average, towards the 1/4 and 
3/4 positions. In the mean time, ter drifts slowly from 
the cell pole towards the midcell, crossing one of the ori 
trajectories. 



E. Application to plasmid segregation — evolutionary 
necessity for segregation strategies. 

Another important biological question concerns plas- 
mid segregation in bacteria. Importantly, most low-copy 
number plasmids, such as Rl, are believed to rely on more 
biological, protein-based active segregation mechanisms 



[e.g., par system ( Garner et al. 2007 Schumacher et al. 
2007)]. Our phase diagram can explain the need for ded- 



icated segregation mechanisms as follows: Since the size 
of typical plasmids is ~10kb, they can be considered as 



Note that, in principle, there are four determining parameters of 
our phase diagram, where the fourth parameter is the number of 
chains n. If we release the constraint imposed on n [e.g., n = 2 
(chromosomes) and n = 4 (chromosomes + plasmids)], while 
keeping constant D and £bulki the two cells are connected by a 
diagonal line in the phase diagram because, then, Rp is the only 
additional variable (Fig. |5| . 

There is an ongoing debate about how long the duplicated ori 
and other loci stay together before splitting during chromosome 
segregation ("cohesion"). See, for example, |Bates and Kleckner 
|| 2005t ; |Na"nninga et q^2002| |; |Nielsen et aZ.| i |2006| ); |Sunako et al 
(20011. 
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FIG. 6 Chromosome segregation in E. coli: simulation vs. experiment. (Left) A series of typical conformations of a replicating 
circular chain from 0% to 99% replication (un-replicated strand in gray, replicated strands in red and blue). We also present 
two sets of segregation pathways (ori-ter trajectories during replication) with and without "replication factory"; the third set of 
simulation where ori starts at the midcell can be found in Supporting Information in Jun et al. | ( 2007 1 . We simulated replication 
factories by enforcing physical proximity of the two replisomes during replication, but we did not fix their position within the 
cell. The dotted lines show the results of 10 individual simulations, and the thick solid lines show the average trajectories 
of ori (red and blue) and ter (black). (Center-to- Right) We juxtapose the simulations with the published data in Bates and 
|Kleckner] ( |2005[ > in an attempt to capture the main features of the experimental observations. For comparison, we used the 
fraction-replicated / as our "universal clock" and scaled the height of the simulation trajectories to match < / < 1 range of 
the data. [Reprinted from Fig. 5 of Jun and Mulder (2006). Copyright (2006) National Academy of Sciences, U.S.A.] 



In C. crescentus, loci trajectories seem even more strik- 
ing (Viollier et al. 2004): one ori stays at the cell 



pole and the duplicated ori moves much faster than the 
growth rate of the cell (Fig. [7]). Indeed, this large dif- 
ference between the cell growth rate and the speed of 
ori is the evidence against Jacob and colleagues' model 
mentioned earlier, at least, in C. crescentus. The rest of 
the loci follow similar trajectories, creating a mirror-like 
image of the organization of duplicated circular chromo- 
somes. Ori is found at the old cell pole and ter at the 
new cell pole; importantly, there is a one-to-one corre- 
spondence between the clock position of the chromoso- 
mal loci and their average physical locations along the 
long axis of the cell. In E. coli, recent data from several 



labs [e.g., ( |Nielsen et al] |2006| [Wang et al] |2006[ )] also 
suggest similar principal linear organization, except that 
their clock positions are rotated by 90 degrees with the 
two chromosome arms between ori and ter occupying 
each cell half. 

To explain these experimental data, we have previously 



proposed the concentric- shell model (Jun and Mulder 



2006), which was inspired by the observation of nucleoid 
compaction in E. coli (see, also, Fig. |8j. Here, our pre- 
diction was that the newly synthesized DNA will move 



much faster in the periphery of the nucleoid near the cell- 
wall membrane than inside the nucleoid body (which is a 
meshwork of chromosomal DNA) , faster than the typical 
timescale of the cell cycle of E. coli (> 20 min.) or C. 
crescentus [» 240 min. in Viollier et al. ( 2004 1] . Indeed, 



this minimal model and assumption reproduced most of 
the major features of the experiment al data in b o th or 
"' For E. coli, we note that Fan et al. 



gamsms. 



(2007) 



to ex- 



have also used a free energy-driven string mode 
plain the experimental data. The most important result 
in their study is that the size of the chromosomal domain 
is important (the larger the domain is, the better chro- 
mosomes segregation), which we have explained above 
using our phase diagram. 

We emphasize that the concentric-shell model is con- 
sistent with other models, as long as they also imply 
the preferential occupation of the nucleoid periphery vol- 
ume by newly synthesized DNA in the early stage of 
cell cycle. These models include the transertion model 



In our previous study ( |jun and~ Muldcr 2006|, this gap was even 
smaller than the width of the chain in our simulations. 
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E 

3= 



Fraction-replicated f 
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FIG. 7 Chromosome segregation of C. crescentus, comparing 
simulation (Left) vs. experiment Viollier el aZ.| ( |2004[ ) (Right). 
The simulated trajectories are the average of 26 individual 
simulation runs (or cells); we show the trajectories of nine 
representative loci (including ori and ter) on the right-arc of 
a circular chromosome for the entire duration of replication 
(up to 99.9%), whereas experimental data are only available 
for trajectories up to 50% of replication. For clarity, we show 
the trajectories only from the onset of replication of each lo- 
cus. A full trajectory of ter is shown, however, to emphasize 
its slow drift from the cell pole to the cell center during repli- 
cation, in contrast to the fast, directed diffusion of ori2 in 
the nucleoid periphery. A main difference from the E.coh 
simulation is the additional assumption that we kept oril in 
the volume near the stalked pole until 10%-20% of the chain 
has been replicated. The concentric-shell model used in the 
simulation is formal and is not inconsistent with other addi- 
tional mechanisms that may act on the directed movement of 
ori2, although in our simulation we did not need any such as- 
sumptions. [Reprinted from Fig. 6 of Jun and Mulder (2006 1. 
Copyright (2006) National Academy of Sciences, U.S.A.] 



by Woldringh (20021, which assumes an interaction be- 



tween the nucleoid and the inner-cell wall membrane via 
co-transcriptional translation and protein translocation, 
or even the hypothetical role of Par proteins on ori trans- 
portation in the cytoplasmic space, namely, outside the 
nucleoid volume (Viollier et al. 2004| ). We believe time- 



lapse experiments with higher spatiotcmporal resolution 
will reveal the nature of directed motion between diffu- 
sive, biased random walk vs. transportation by polymer- 
ization or motor proteins. 



III. DISCUSSION AND FREQUENTLY ASKED 
QUESTIONS 




Topologically indepen- 
dent structural unit 



close-packing in a cell 




FIG. 8 Illustration of chromosome as a close packing of a 
string of Pincus blobs (see, also, Appendix. [b| . The blobs 
schematically represent correlation-length structural units of 
the chromosome, which can be formed by stretching (and 
twisting) of the chain. (Although not essential to our the- 
ory presented in this article, independence of the topologi- 
cal states of individual blobs may be achieved by nucleoid- 
associated proteins.) To maintain an average linear organi- 
zation, we envision the blobs forming a helical coil that fills 
the nucleoid volume. If the helicity reverses, nucleoid enve- 
lope may develop constricted, bi-lobed or multi-lobed shape. 
Molecular crowding and possible interactions between the in- 
ner cell-wall membrane and dsDNA protruding from the nu- 
cleoid is not shown for clarity. The gap between the nucleoid 
and the inner cell-wall membrane is not drawn to scale and 
may vary from organism to organism. 



Based on the piston model of the bacterial cell, we have 
presented the above physical biology model of bacterial 
chromosomes. Importantly, our model makes an exper- 
imentally testable predictions whether duplicated chro- 
mosomes will mix or segregate, and we have critically 
examined our predictions against the published data of 
E. coli B/r strain. 

We emphasize the importance of the concept of corre- 
lation length £, a central parameter of our chromosome 
model. To recapitulate, this is the size of the chromoso- 



mal structural unit within which the density fluctuations 
(and, loosely, the motions) of the chain segments are cor- 
related. The major assumption underlying our model is 
that, although a real bacterial chromosome is a very com- 
plex entity because of supercoiling and the presence of 
various nucleoid-associated proteins (whose roles are still 
largely unknown), there is a single length scale £ that 
characterizes the nucleoid. Then, the chromosome can 
be interpreted as an entity consisting of these structural 
units, and chromosome organization in bacteria can be 
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understood as a close packing of a string of Pincus blobs 
We illustrate this view in Fig. [8] and leave the scaling 
analysis for equivalence between the Pincus chain (Pin- 
1976 1 and the compressed chain in Appendix. [B] 



cus 



Based on our results discussed in the previous sections, 
there are two general principles in understanding chromo- 
some organization and segregation in bacteria. 

I. The larger the correlation length is, the better 
chromosomes segregate. 

II. The larger the aspect ratio is (i.e., the longer 
and/or the narrower the cell is), the more 
chromosomes tend to demix. 



Cell growth and division. Since the cell size (relative to 
the correlation length) is also important, proteins that 
regulate cell growth and division are also important to 
ensure that the cell reaches the appropriate size or mass 
for proper segregation |Weart et al. 2007[ ). In this con- 
text, the length-measuring "devices" such as MinCDE 
are also important, and proteins involved in cell division 
such as FtsZ are also relevant [see ( Lutkenhaus[ 



and references therein], if cell constriction should help 



2007) 



resolve partially overlapping nucleoids ( Huls et al. 1999 1 



B. Shouldn't a hypothetical motor protein be enough to 
segregate chromosomes in bacteria? 



Below, we discuss various questions concerning our 
chromosome model. 



A. Role of proteins 

Proteins directly and indirectly change the physical 
properties of chromosomes and the cell morphology and, 
thus, they are important for chromosome organization 
and segregation in the following three contexts. 

Correlation length of chromosome. Various SMC or 
nucleoid-associated proteins such as MukBEF, HU, FI- 
NS, IHF and gyrase may change the size of the chromo- 
somal structural unit and, thus, can help ensure success- 



ful chromosome segregation (Stavans and Oppenheim 
[20061). From our point of view, the best evidence has been 



given by 


Sawitzke and Austin 


(2000 


1 [see, also, 


Holmes 


and Cozzarclli (2000|) and 


Dasgupta et al. 


(20001 


. These 



authors have demonstrated that severity of disruption 
of chromosome segregation of muk~ E.coli can be con- 
trolled by the level of supercoiling in the cell. Thus, 
neither class of proteins is the dedicated segregation ma- 
chinery, but their role can be understood as, within the 
theoretical framework presented here, to increase the ef- 
fective correlation length of the chromosome above the 
critical size (Eq. |3| , which is sufficient for proper segre- 
gation (see I above). 

We also note that, because of supercoiling, a more re- 
alistic coarse-grained topology of a bacterial chromosome 
is that of the branched polymer rather than a purely 
linear chain. This branched structure will enhance the 



tendency of chromosome pushing even more (Jun and 
Mulderl 120061 iTrun and Markol 119981 IVilgisl 12000ft. 



Cell shape, MreB, and bacterial 11 cyto skeletal" proteins. 



MreB and other proteins ( Gerdes et al. 2004 1 are 
important to maintain the high aspect ratio of the 
cell (see II above). It is our contention that, if these 
proteins contribute to chromosome segregation, it is 
about changes in cell shape as our model suggests. 



No. There is increasing experimental evidence that 
ori and perhaps some other chromosomal loci are local- 
ized or te thered at specific intracellular po sitions [see, for 
example, Stavans and Oppenheim (2006) and references 
therein]. However, one still must explain how the hy- 
pothetical transportation of a small fraction (e.g., ori) 
of millions of basepairs of DNA from one position inside 
the cell to another can dictate, if any at all, the directed 
movement, segregation and organization of the rest of 
the chromosome. These observations thus require under- 
standing of more basic physical principles (e.g., Fig. [4]). 

We note that there are special cases where motor 
proteins are indeed involved in transportation of DNA 
from one position to another inside a bacterial cell. The 
examples include SpoIIIE for sporulation in B.subtilis 
and FtsK for dimcr resolution and other "rescue" tasks 
in E.coli Barre (20071 and references therein]. However, 



it is important to realize that these proteins translocate 
DNA, i.e., they take advantage of the directionality 
provided by the septum, which is entirely different from 
segregation of replicating chromosomes. 



C. What is the timescale? 

In the segregation regime, segregation is a driven pro- 
cess and two int ermingled chains drift at typical time 



scal e of t ~ N 2 ( |Arnold and Jun[ |2007| | Jun and MuT 



der 



2006J, which is much faster than diffusion (repta- 
tion) timescale of r ~ N 3, in the channel. Timescale 
of polymer motions in the presence of confinement in- 

and 



2007) 



volves severe finite-size effects (Arnold et al. 
it is far more difficult to predict the timescale for seg- 
regation of duplicating chains. For the concentric-shell 
model in our simulations for E.coli and C. crescentus, 
we assumed a separation of timescale that the relaxation 
time of the replicating segment of DNA is much smaller 
than the typical timescale of cell cycle (Sec. II.Fl. Such 
fast movements of DNA segment have indeed been re- 
ported in recent experiments, that the 15kb minimalized 
plasmid RK2 explores about 50% of the cell volume only 



within 2 minutes ( |Derman et al. 2008 ) 
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(a) 



D 



(b) 



Khokhlov[ |1994| |Jun et al] |2007| |Sakaue and Raphael 




H D' 




FIG. 9 Two possible scenarios of the role of molecular crowd- 
ing on chain interactions, (a) segregation (b) mixing. Each 
chain consists of N monomers. Calculations show that crowd- 
ing, in the first approximation, should not influence mixing 
vs. de-mixing of the chains (see text). 



We do feel that the issue of timescale can be addressed 
adequately only via experiments using novel physical 
techniques. 



D. How about topological^ catenated chromosomes? 

Imagine two topologically concatenated ring polymers 
confined in a rectangular box, where our phase diagram 
in Fig. [4] predicts segregation of two linear chains of 
the same contour length as the rings. These two ring 
polymers will still segregate to occupy each half of the 
box, i.e., the topology of the system localizes on average 
at the center of the box (unpublished results). It is an 
intriguing question whether this localization of chain 
topology, which is a global property of the system, 
will influence the action of topoisomerases, which have 
access to only local information. We speculate that 
the entropy-driven directionality of the movements and 
relative orientation of the confined chains will, at least 
in part, help topoisomerases decatenate the replicated 
chains before cell division. 



E. Bacterial cells are very crowded: Will molecular 
crowding influence chromosome organization and 
segregation? 

In the first approximation, no. One may foresee two 
opposing arguments: (A) As the chains are compressed 
by depletion effect, the inner space occupied by the chains 
has higher density of polymers, i.e., less accessible vol- 
ume by the other chain and, thus, the tendency of de- 
mixing will increase because of squeezing, or (B) Molec- 
ular crowding exerts an effective osmotic pressure, which 
will make the chains mix. 

In fact, a simple calculation suggests that there is a 
balance between these two effects. To see this, let us con- 
sider the free energy of the compressed chains in Fig. [9 



as follows (Cacciuto and Luijten 2006 Grosberg and 



2006) 



/3F A = 2C 



13Tb 



C 



(2NY 
D' 



(4) 
(5) 



where C is the pref actor. These two free energies for 
the compressed chains become equal at D' = 2 1//3 D, 
i.e., when their total volumes are the same, Va = Vb- 
Thus, compression due to depletion by itself will not in- 
fluence mixing vs. de-mixing of the chains, even if it 
could bring the chains together and change their enve- 
lope shape (surface tension). 10 Computer simulation re- 
sults also support the neutral effect of molecular crowding 
(Axel Arnold, personal communication). 

On the other hand, crowding can influence the local 
organization of chromosome. For instance, looping is 
one way to achieve compaction of chromosomes, and the 
entropy gain by depletion attraction between DNA seg- 
ments can be larger than the entropy loss by DNA loop- 
ing (Marenduzzo et al. 2006). 



F. Bacteria exist in various cell shapes and composition of 
chromosomes. Are there general strategies for successful 
chromosome segregation in bacteria? 

Since chromosome segregation is one of the defining 
processes of any cell, its basic mechanism must have 
been conserved across branches of life and organisms of 
diverse shapes. Based on our phase diagram (Fig. [4} 
and the two general rules I Sz II presented above, we 
can speculate about how bacteria may create favorable 
physical conditions for partitioning the chromosomes 



(Fig. 101 



Filamentous bacteria. We propose an entropy- 
driven, "random" segregation mechanism for filamentous 
bacteria. As we have shown both analytically and nu- 
merically, polymers confined in a narrow cylindrical ge- 
ometry strongly resist overlap and, thus, repel one an- 



other (Arnold and Jun 2007 Jun and Mulder 2006), 



where the timescale of disentanglement of overlapping 
polymers is much shorter than that of pure diffusion. 
Indeed, recent experimental study on the filamentous 



10 The origi n of this result i s due t o the form of the free energy 
in Eq. |4| ||Grosberg et al.\ |1982| |Jun et al.\ |2007| |Sakaue and] 
Raphael) |2006) . which has been constructed to be a function 
of only the monomer density. In other words, regardless of 
the shape of the envelope surrounding the chains, the free en- 
ergies will be the same as long as the volumes are also the same, 
and vice versa. Additional consideration of the interactions at 
the surface will change the envelope shape, but not segrega- 
tion/mixing. 
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(a) 



□ © 



(b) symmetry breaking of cell shape 




square bacteria ("2-dimension") 



spherical bacteria 



FIG. 10 Possible cell-shape-dependent strategies for chromo- 
some segregation, (a) According to our phase diagram, fil- 
amentous and rod-shaped cells are most favored by entropy- 
driven segregation, (b) Thin, square bacteria such as Walsby's 
square bacterium change their direction of growth after every 
division (typically they are found in composition of lxl, 2x2, 
4x4, 8x8 cells and so on). On the other hand, round cells 
divide at alternating perpendicular planes. The arrows in the 
round cells indicate oscillation of the MinCDE proteins along 
the long axis of the ellipsoidal cells, which help determine the 
mid plane of the cell. 



cyanobacterium Anabaena sp. PCC 7120 has revealed 
Gaussian distribution of DNA content in each daugh- 
ter cells after the septum formation (IHu et al. 



2007). 



Importantly, the variation of the distribution was much 
larger than a value expected if the two daughter cells were 
identical, suggesting a random event involved in chromo- 
some segregation. Also, Anabaena sp. PCC 7120, like 
many other filamentous bacteria, is polyploid and con- 
tains multicopy of chromosomes per cell, ~ 10 (Hu et al. 



2007 ), which seems to moderate the effect of random seg- 



regation. Since chromosomes occupy much larger volume 
than plasmids inside the cell, this apparent random seg- 
regation process is sufficient for the viability of the cell 
as long as multiple copies of chromosomes are produced 
before division. 

Rod-shaped bacteria. We have critically examined the 
model system E. coli B/r strain and explained why 
duplicated chromosomes will not mix. Thus, other 
organisms of similar cell shape and volume must benefit 
from the entropy-driven driving force of segregation, 
with or without any organism-specific segregation mech- 
anisms, as long as the the correlation of chromosome is 
also comparable to that of E. coli (« 100 nm). 

Square bacteria. Perhaps one of the most striking bacte- 



rial shapes is that of Walsby's square bacterium ( jWalsby 
2000 ) . This organism has a very thin and square-shaped 
cell (sa lOOnm thick and ~ 2[im wide with aspect ratio 
1), a "2-dimensional" creature resembling a postage 
stamp. As we have explained above, polymers in highly 
symmetric confinement can readily intermingle. In other 
words, entropy- favored condition for segregation can be 
created when the cell breaks its symmetry during growth 
or division. Indeed, microscopy images so far suggest 



that that square bacteria grow and divide at alternating 
perpendicular axes and planes, respectively [Fig. |10[ b)]. 
Also, in a thin slab (i.e., 2D confinement), using the 
approach taken in I Jun et al. p007|, it is straightforward 



to show that polymers repel much more strongly than 
in bulk 3D. 11 (Note that the cells are typically observed 
in 2x2, 4x4, 8x8 stamp-like configurations.) Therefore, 
we believe that simple symmetry breaking by cell 
growth and the thin 2D geometry suffice to separate the 
chromosomes entropically. 

Spherical bacteria. The perfect symmetry of the cell 
shape means that the confined chains do not have any 
preferred conformations between mixing and de-mixing, 
although their global reorganization could readily be 
achieved (Jun et al. 20071. Since little data is available 



on chromosome organization in spherical bacteria, we 
can only speculate about possible contributing factors 
to segregation. At the chromosome level, supcrcoiling- 
induced branched structures of the chromosome will 



increase the tendency of demixing (Jun and Mulder 



2006 Trun and Marko 1998 Vilgis 20001, perhaps 



helped by nucleoid-associated proteins to increase the 
correlation length (rule I above). More importantly, 
at the cell level, symmetry breaking of the cell shape 
and invagination of the cell during division could help 
resolve partially intermingled chromosomes [Fig. [Toj^b)] . 
Indeed, real cells are never perfectly spherical, and 



Huang and Wingreenj (20041 have shown numerically 



that Min-protein oscillations can be achieved along the 
long axis of the nearly round cell, where the equatorial 
radii differ by as small as 5%. This may explain the 
observed division of round cells at alternating perpen- 
dicular planes (Corbin et al. 20021. Along with the 



case of the filamentous bacteria discussed above, it is 
tempting to speculate that chromosome segregation in 
spherical bacteria is also a "random" process and, thus, 
polyploidy can be a strategy to increase the chance 
of successful distribution of chromosomes to daughter 
cells. For example, the endosymbiotic bacteria of 
aphids, Buchnera, which cannot divide by itself outside 
the host eukaryotic cell, contains over 100 copies of 
genome (Komaki and Ishikawa 19991. 



G. Early life and artificial cells 

Previously, we proposed that the entropy-driven seg- 
regation of polymers could be implemented in designing 
an artificial cell, and that it has implications on early 



life (Jun and Mulder 2006). How could an artificial cell 



11 The free-energy cost of chain overlapping crosses over from /3 T ~ 
n dv d /(du d -l) to ^ n d/(d-2) j w h ere u d = 3/{d + 2 ) is th e Flory 
exponent and d the spatial dimension. [Jun et aT] | |2007| l; B.-Y. 
Ha, personal communication] 
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Division 




replication and/or 
polymerization 



Polymer repulsion S^B* 1 
elongation of membrane 

FIG. 11 Entropy-driven segregation of encapsulated polymers 
and its potential role on the cell cycle of a protocell. Repul- 
sion between the replicated chains can break the spherical 
symmetry of the vesicle, and, thus, may facilitate its spon- 
taneous division. Moreover, the daughter vesicles containing 
the duplicated polymers may have similar sizes. (The three 
minisketches at the top represent input from environment for 
feeding the growth.) These scenarios could be tested both 
theoretically and experimentally. [Adapted from Figure 3 in 
|Szostak et al\ ( [200T| .] 



quantitative sciences are making significant contribu- 
tions. From a physicist's point of view, however, it is 
interesting and important to notice that many of the con- 
tributions made by physicists so far have involved bring- 
ing new tools and ways of thinking to (systems) biology, 
rather than new understanding about the role that ba- 
sic physical principles play in governing the fundamental 
biological processes inside the cell. But aren't biological 
entities also physical objects? The membranes defining 
the envelope of the organelles; the polymers constituting 
proteins, spindles, and chromosomes - these are physical 
objects and, although biophysical approach has uncov- 
ered many of their important physical properties, their 
biological implications need to be more fully explored. 

Thus, in seeking a more direct relationship between 
physics and biology, with the example of bacterial chro- 
mosome segregation, I hope to have convinced the reader 
that there is a seamless "dialogue" between physical and 
biological processes involved in the bacterial cell cycle; 
and that deeper understanding of the relationship be- 
tween the physical principles and their biological impli- 
cations may even shed new light on the major tran sitions 
in evolution ( Maynard-Smith and Szathmary 1998). 12 
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achieve its basic cellular processes? While Szostak and 
colleagues have proposed a system of a self-replicatin g 
vesicle and replicases as a protocell dSzostak et al. 2001 1, 



cell division mechanism as well as the size regulation of 
daughter cells of a protocell are far from being under- 
stood. 

We suggest that two long polymers, formed by replica- 
tion/ligation of smaller molecules within a spherical vesi- 
cle repel each other because of entropy. If the polymer- 
polymer repulsion is strong enough to break the spher- 
ical symmetry of the vesicle, this process may lead to 
membrane fission, preferentially at the mid-cell position, 
which is defined by the surface of contact between the two 
repelling polymers. This would further help regulate the 
size distribution of newly formed daughter cell as illus- 
trated in Fig. 11 Although it is in general a formidable 



theoretical problem to predict the shape of the closed 
membrane in a given environment, it would be highly 
desirable to obtain a phase diagram similar to Fig. [4] for 
soft walls, and possibly test our proposed idea here exper- 
imentally using macromolecules encapsulated in vesicles 
or micro drops. 



IV. FINAL REMARKS 

In recent years, biology has become increasingly inter- 
disciplinary. In particular, in the new discipline of sys- 
tems biology, researchers trained in physical and other 
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12 For instance, in the presence of multiple chromosomes in cu- 
karyotes, the cell faces an additional challenge in chromosome 
segregation because each daughter cell must receive not only the 
correct number but also the correct set of chromosomes. From 
out point of view, this strongly argues why mitosis requires more 
sophisticated segregation machinery as well as the checkpoints; 
whereas entropy might be sufficient for chromosome segregation 
in most bacteria. 

On the other hand, the hallmark of mitosis is the compact 
sister chromosomes of well-defined shapes being held together, 
waiting for the mitotic spindle to separate them. However, these 
chromosomes must go through a period of intermingling during 
replication, before mitosis. Thus, there is a process analogous to 
bacterial chromosome segregation, namely, de-mixing of repli- 
cated eukaryotic sister chromatids, which we distinguish from 
segregation by spindle. Although the exact mechanism of this 
de-mixing process has remained elusive, we believe that the basic 
physics behind is the same as what we presented in this article. 
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APPENDIX A: More technical description of the phase 
diagram in Fig. [4j 

Below, we provide more technical descriptions of indi- 
vidual regimes with references so that the reader can 
reproduce the boundary conditions of the phase diagram 



where V = L/R < 2, R' = R/Rq < 1, D' = D/R Q . 
Since dpT tot /dR' = Zk/2 > at R' = L'/2 (namely, 

indeed 



when each chain occupies each cell half), Eq. A2 
predicts a free-energy barrier against mixing of two 
segregated chains. Moreover, one can show numerically 



in Fig. |4| [See, also, 


Brochard and de Gennes| 


1979 


I; function of R 1 > L'/2 for V > 0.8. 13 
]■ 


Daoud and de Gennes 


(1977 


); 


Teraoka and Wang|( 


2004 



I. This is the "bulk" dilute regime, where the interchain 
distances are much larger than the size of the chain in 
confinement. When D = L = Rp . the correlation length 
£buik a l so equals the size of the chain Rp, and the two 
chains can overlap at the free-energy cost of order ksT 
[Fig.ga)] ( |Grosberg efalj [l982| [Tun et al. [[2007] . 



II. Imagine dilute solution in a long cylindrical box, 
where there is a free space between the two chains 
[Fig. Bib)]. The size of individual chains scales as i?o ~ 
(N/ g)D ~ ND~ 2 / 3 (where g is the number of monomers 
per blob). In case of overlapping, the chains strongly 
repel one another and, thus, their segregation is much 
faster than reptation ( typical timescale of r 



N 2 



vs. 



r ~ A^ 3 , respectively) (Arnold and Jun 



2007 



Jun and 



Mulder 20061. Note that £buik > D in this regime, 
whereas the blob size of the chain due to confinement 
is D. Thus, the monomer- monomer correlation length 
inside the box is £box — D < £buik- 

If we gradually reduce the aspect ratio of the box, 
keeping constant the volume of the box (thus, fixing 
£tmik) 5 the number of blobs of each chain decreases and 
eventually the two chains make contact at L — 2Rq, i.e., 
Cbuik = £box = D (y = x). This situation is depicted 
in Fig. |4jc), at which the chains enter the semi-dilute 
regime inside the elongated box. 

III. If we continue the above process of reducing the as- 
pect ratio of the box, the two chains are gradually com- 
pressed but do not mix as long as their principal con- 
formations are linear [Fig. |4jd).] As we discussed in 
Sec. II. D this is the most important regime for rod- 
shaped bacteria. To see the segregation of overlapping 
chains, let us employ the following Flory-type free en- 
ergy of individual chains: 



R 2 



(N/g)D* 



D(N/g) 2 
R 



(Al) 



IV, V, VI (mixing regimes). As the aspect ratio ap- 
proaches 1, the two chains lose their linear conformations 
and mix with each other. The reason for this transition 
is that the bulk correlation length £buik becomes much 
smaller than both D and L, and the chain conformation 
becomes a random walk of a string of blobs (blob size 
£b uik)- We d e scribe d the physics of this regime in detail 
in 



Jun et al. (20071. 



At the boundary between III and IV, transition oc- 
curs because individual chains are contracted and lose 
their linear ordering. The size of individual chain is 
that of the random- walk conformation, R» ~ jV 1 / 2 ^ 1 / 8 
with the volume fraction of the monomer <f> ~ N/R^D. 
Based on t his, it is easy to show the boundary condition 

y 



regime is given by 



12//7 (Lai et al. 1997). Since the free energy in this 



IV 



D 2 L 

^bulk 



(A3) 



(Cacciuto and Luijten 


2006 Grosberg and Khokhlov 


1994 


Jun et al. 


2007 


Sakaue and Raphael 


2006 


1, the 



above boundary condition is translated to D = k ■ £buik, 
which is in symmetry with the condition L = k ■ D. Sim- 
ilarly, it is straightforward to show that the boundary 
condition k = 1 between IV and V is identical to y — x 9 ^ 4 
[Fig- He)]. 

In regime V, the shape of the box resembles a thick 
slab. As the thickness of the slab reaches the size of the 
blob, i.e., £buik = L = k ■ D [Fig. |4je)], the two chains 
enter the dilute regime in a thin, closed slab [Fig. |4jg)], 
which is the dual regime of II. Thus, the transition be- 
tween V and VI occurs when the chains start to feel the 
size of the slab, i.e., at D ~ Rn — iV 1 / 2 ^ -1 / 8 , where the 



volume fraction is given by $ ~ ff/^bulk' w ith g 



where (3 = l/k B T ( jJun et a?.||2008[ ). On the other hand, We thus obtain y = x 6 (L = k ■ D) 



Al" 

?bulk- 



it costs ~ fc^T of free energy for two blobs to significantly 



overlap regardless of their size ( 


Grosberg et al. 


1982| Jun 


et al. 2007 


Jun and Mulder 2006|). Thus, the total free 



energy of the system of two partially overlapping chains 
described in Fig. |4jh) can be written as 



R 



0^1 OR) 





' 4 




B' 2 -U 1 



















(A2) 



13 Since the N 2 /R dependence in Eq. |Al| is an well-known overes- 
timate of Flory-type approach for R <g _Rqi the L' pb 0.8 is an 
upperbound. In other words, two chains should segregate in an 
even shorter box. For the same reason, we cannot use Eq. |A2| 
to determine the condition for crossover to regime IV. To our 
knowledge, a correct form of the free energy that is valid even 
for largely perturbed chains (namely, for all range of R) is not 
available. 
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APPENDIX B: Self-consistent mapping between the Pincus 
and the close-packed chains. 

In the Pincus picture of a stretched chain, the chain 
breaks up into a series of blobs of size £ (Fig. [sj> . The free 
energy cost can be estimated by counting the number of 
blobs, 



N 
9 



R 



N 



(Bl) 



where g ~ t^-l" is the number of monomers per blob and 
R the end-to-end distance of the chain. 

Let us now imagine close-packing of the chain consist- 
ing of blobs of size £ in a cylinder of diameter P. If this 
process does not cost any additional free energy, using 



Eq. Bl we can self-consistently obtain the confinement 
free energy accordingly as follows. 



LP 2 



N 



LP 2 



(B2) 



where L is the longitudinal size of the close-packed chain 



in the cylinder. Eq. B2 has the same scaling form of the 
free energy cost for confining a self-avoiding chain in a 
volume V , 



V 

#3 



(B3) 



(see Eq. A3 1 , and we thus have a self-consistent picture 



for interpreting compressed chain based on the Pincus 
chain. 

The essence of our analysis here is that the free energy 
for deformation costs ~ fc^T/blob regardless of the na- 
ture of deformation and, thus, is mappable one another 
(e.g., between stretching vs. confinement). 



References 

Adachi, S., K. Hori, and S. Hiraga, 2006, "Subcellular Posi- 
tioning of F Plasmid Mediated by Dynamic Localization of 
SopA and SopB," J Mol Biol 356, 850-863. 

Alder, B. J., and T. E. Wainwright, 1957, "Phase Transition 
for a Hard Sphere System," J Chem Phys 27, 12081209. 

Arnold, A., B. Borzorgui, D. Frenkel, B.-Y. Ha, and S. Jun, 
2007, "Unexpected relaxation dynamics of a self-avoiding 
polymer in cylindrical confinement," J Chem Phys 127, 
164903. 

Arnold, A., and S. Jun, 2007, "Timescale of entropic seg- 
regation of flexible polymers in confinement: Implications 
for chromosome segregation in filamentous bacteria," Phys 
Rev E 76, 031901. 

Barre, F.-X., 2007, "FtsK and SpoIIIE: the tale of the con- 
served tails," Mol Microbiol 66, 1051-1055. 

Bates, D., and N. Kleckner, 2005, "Chromosome and Repli- 
some Dynamics in E. coli: Loss of Sister Cohesion Triggers 
Global Chromosome Movement and Mediates Chromosome 
Segregation," Cell 121, 899-911. 



Berkmen, M. B., and A. D. Grossman, 2006, "Spatial and 

temporal organization of the Bacillus subtilis replication 

cycle," Mol Microbiol 62, 57-71. 
Brenner, M., and J. -I. Tomizawa, 1991, "Quantitation of 

ColEl-encoded replication elements," Proc. Nat. Acad. Sci. 

88, 405-409. 

Brochard, F., and P. G. de Gennes, 1979, "Conformation of 
molten polymers inside small pores," J. Phys. (Paris) 40, 
L399-L401. 

Cacciuto, A., and E. Luijten, 2006, "Self- Avoiding Flexible 
Polymers under Spherical Confinement," Nano Lett 6, 901- 
905. 

Corbin, B. D., X.-C. Yu, and W. Margolin, 2002, "Exploring 
intracellular space: function of the Min system in round- 
shaped Escherichia coli," EMBO J. 21, 1998-2008. 

Cunha, S., T. Odijk, E. Siileymanoglua, and C. L. 
Woldringh, 2001, "Isolation of the Escherichia coli nu- 
cleoid," Biochimie 83, 149-154. 

Danchin, A., P. Guerdoux-Jamet, I. Moszer, and P. Nitschke, 
2000, "Mapping the bacterial cell architecture into the 
chromosome," Phil. Trans. R. Soc, Lond. B 355, 179-190. 

Daoud, M., and P. G. de Gennes, 1977, "Statistics of macro- 
molecular solutions trapped in small pores," J. Phys. 
(Paris) 38, 85-93. 

Dasgupta, S., S. Maisnier-Patin, and K. Nordstrom, 2000, 
"New genes with old modus operandi: The connection be- 
tween supercoiling and partitioning of DNA in Escherichia 
coli," EMBO Reports 1, 323-327. 

Deng, S., R. A. Stein, and N. P. Higgins*, 2005, "Organization 
of supercoil domains and their reorganization by transcrip- 
tion," Mol Microbiol 57, 1511-1521. 

Derman, A. I., G. Lim-Fong, and J. Pogliano, 2008, "Intracel- 
lular mobility of plasmid DNA is limited by the ParA family 
of partitioning systems," Mol Microbiol 67, 935-946. 

Doi, M., and S. F. Edwards, 1986, The Theory of Polymer 
Dynamics (Oxford University Press, Oxford, UK). 

Dworkin, J., and R. Losick, 2002, "Does RNA polymerase 
help drive chromosome segregation in bacteria?," Proc Nat 
Acad Sci 99, 14089-14094. 

Elmore, S., M. Miiller, N. Vischer, T. Odijk, and C. L. 
Woldringh, 2005, "Single-particle tracking of oroC-GFP 
fluorescent spots during chromosome segregation in Es- 
cherichia coli," J Struct Biol 151, 275-287. 

Espeli, O., and F. Boccard, 2006, "Organization of the Es- 
cherichia coli chromosome into macrodomains and its pos- 
sible functional implications," J Struct Biol 156, 304-310. 

Fan, J., K. Tuncay, and P. J. Ortoleva, 2007, "Chromosome 
segregation in Escherichia coli division: A free energy- 
driven string model," Comp Biol Chem 31, 257-264. 

Fiebig, A., K. Keren, and J. A. Theriot, 2006, "Fine-scale 
time-lapse analysis of the biphasic, dynamic behaviour of 
the two Vibrio cholerae chromosomes," Mol Microbiol 60, 
1164-1178. 

Flory, P., 1953, Principles of Polymer Chemistry (Cornell 

University Press). 
Fogel, M. A., and M. K. Waldor, 2006, "A dynamic, mitotic- 

like mechanism for bacterial chromosome segregation," 

Genes Dev. 20, 3269-3282. 
Garner, E. C, C. S. Campbell, D. B. Weibel, and R. D. 

Mullins, 2007, "Reconstitution of DNA Segregation Driven 

by Assembly of a Prokaryotic Actin Homolog," Science 

315, 1270-1274. 
de Gennes, P.-G., 1979, Scaling Concepts in Polymer Physics 

(Cornell University Press, Ithaca, NY). 



15 



Gerdes, K., J. Moller- Jensen, G. Ebersbach, T. Kruse, and 
K. Nordstrom, 2004, "Bacterial Mitotic Machineries," Cell 
116, 359-366. 

Grosberg, A. Y., P. G. Khalatur, and A. R. Khokhloo, 
1982, "Polymeric Coils with Excluded Volume in Dilute 
Solution: the Invalidity of the Model of Impenetrable 
Spheres and the Influence of Excluded Volume on the 
Rates of Diffusion-Controlled Intermacromolecular Reac- 
tions," Makromol. Chem., Rapid Commun. 3, 709-713. 

Grosberg, A. Y., and A. R. Khokhlov, 1994, Statistical 
Physics of Macromolecules (AIP Press) . 

Holmes, V. F., and N. R. Cozzarelli, 2000, "Closing the ring: 
Links between SMC proteins and chromosome partitioning, 
condensation, and supercoiling," Proc Nat Acad Sci 97, 
13221324. 

Hu, B., G. Yang, W. Zhao, Y. Zhang, and J. Zhao, 2007, 
"MreB is important for cell shape but not for chromosome 
segregation of the filamentous cyanobacterium Anabaena 
sp. PCC 7120," Mol Microbiol 63, 1640-1652. 

Huang, K. C, and N. S. Wingreen, 2004, "Min-protein oscil- 
lations in round bacteria," Phys. Biol. 1, 229-235. 

Huls, P. G., N. O. E. Vischer, and C. L. Woldringh, 1999, 
"Delayed nucleoid segregation in Escherichia coli," Mol Mi- 
crobiol 33, 959-970. 

Jacob, F., S. Brenner, and F. Cuzin, 1963, "On the Regulation 
of DNA Replication in Bacteria," J Mol Biol 28, 329-348. 

Jun, S., A. Arnold, and B.-Y. Ha, 2007, "Confined Space and 
Effective Interactions of Multiple Self- Avoiding Chains," 
Phys Rev Lett 98, 128303. 

Jun, S., and B. Mulder, 2006, "Entropy-driven spatial orga- 
nization of highly confined polymers: Lessons for the bac- 
terial chromosome," Proc Nat Acad Sci 103, 12388-12393. 

Jun, S., D. Thirumalai, and B.-Y. Ha, 2008, "Compression 
and stretching of a self-avoiding chain under cylindrical 
confinement," Phys Rev Lett (in press) . 

Kleckner, N., D. Zickler, G. H. Jones, J. Dekker, R. Padmore, 
J. Henle, and J. Hutchinson, 2004, "A mechanical basis for 
chromosome function," Proc Nat Acad Sci 101, 12592- 
12597. 

Komaki, K., and H. Ishikawa, 1999, "Intracellular Bacterial 
Symbionts of Aphids Possess Many Genomic Copies per 
Bacterium," J Mol Evol 48, 717-722. 

Kubitschek, H. E., 1981, "Bilinear Cell Growth of Escherichia 
coli," J Bacteriol 148, 730-733. 

Lai, J., S. K. Sinha, and L. Auvray, 1997, "Structure of Poly- 
mer Chains Confined in Vycor," J. Phys. II (France) 7, 
1597-1615. 

Lemon, K. P., and A. D. Grossman, 2001, "The extrusion- 
capture model for chromosome partitioning in bacteria," 
Genes Dev 15, 2031-2041. 

Lutkenhaus, J., 2007, "Assembly Dynamics of the Bacterial 
MinCDE System and Spatial Regulation of the Z Ring," 
Annu Rev Biochem 76, 539-562. 

Marenduzzo, D., C. Micheletti, and P. R. Cook, 2006, 
"Entropy-Driven Genome Organization," Biophys J 90, 
3712-3721. 

Marko, J. F., and E. D. Siggia, 1995, "Statistical mechanics 

of supercoiled DNA," Phys Rev E 52, 2912-2938. 
Marko, J. F., and E. D. Siggia, 1997, "Polymer Models of Mei- 

otic and Mitotic Chromosomes," Mol Biol Cell 8, 22172231. 
Maynard-Smith, J., and E. Szathmary, 1998, The Major 

Transitions in Evolution (Oxford University Press, USA). 
Nanninga, N., M. Roos, and C. L. Woldringh, 2002, "Models 

on stickiness of replicated Escherichia coli oriC," Microbi- 



ology 148, 3327-3328. 

Nielsen, H. J., J. R. Ottesen, B. Youngren, S. J. Austin, and 
F. G. Hansen, 2006, "The Escherichia coli chromosome is 
organized with the left and right chromosome arms in sep- 
arate cell halves," Mol Microbiol 62, 331-338. 

Nordstrom, K., and K. Gerdes, 2003, "Clustering versus ran- 
dom segregation of plasmids lacking a partitioning func- 
tion: a plasmid paradox?," Plasmid 50, 95-101. 

Onsager, L., 1949, "The effect of shape on the interaction of 
colloidal particles," Ann. NY. Acad. Sci. 51, 627659. 

Pincus, P., 1976, "Excluded Volume Effects and Stretched 
Polymer Chains," Macromolecules 9, 386-388. 

Romantsov, T., I. Fishov, and O. Krichevsky, 2007, "In- 
ternal Structure and Dynamics of Isolated Escherichia 
coli Nucleoids Assessed by Fluorescence Correlation Spec- 
troscopy," Biophys J 92, 2875-2884. 

Sakaue, T., and E. Raphael, 2006, "Polymer Chains in 
Confined Spaces and Flow-Injection Problems: Some Re- 
marks," Macromolecules 39, 2621-2628. 

Sawitzke, J. A., and S. Austin, 2000, "Suppression of chromo- 
some segregation defects of Escherichia coli muk mutants 
by mutations in topoisomerase I," Proc Nat Acad Sci 97, 
1671-1676. 

Schrodinger, E., 1944, What is life? (Cambridge University 
Press, Cambridge, UK). 

Schumacher, M. A., T. C. Glover, A. J. Brzoska, S. O. Jensen, 
T. D. Dunham, R. A. Skurray, and N. Firth, 2007, "Segro- 
some structure revealed by a complex of ParR with cen- 
tromere DNA," Nature 450, 1268-1271. 

Stavans, J., and A. Oppenheim, 2006, "DNAprotein interac- 
tions and bacterial chromosome architecture," Phys Biol 3, 
R1-R10. 

Sunako, Y., T. Onogi, and S. Hiraga, 2001, "Sister chro- 
mosome cohesion of Escherichia coli," Mol Microbiol 42, 
1233-1241. 

Szostak, J. W., D. P. Bartel, and P. L. Luisi, 2001, "Synthe- 
sizing life," Nature 409, 387-390. 

Teleman, A. A., P. L. Graumann, D. C.-H. Lin, A. D. Gross- 
man, and R. Losick, 1998, "Chromosome arrangement 
within a bacterium," Curr Biol 8, 1102-1109. 

Teraoka, I., and Y. Wang, 2004, "Computer simulation stud- 
ies on overlapping polymer chains confined in narrow chan- 
nels," Polymer 45, 3835-3843. 

Trun, N. J., and J. F. Marko, 1998, "Architecture of a Bac- 
terial Chromosome," ASM News 64, 276-283. 

Vilgis, T. A., 2000, "POLYMER THEORY: PATH INTE- 
GRALS AND SCALING," Phys Rep 336, 167-254. 

Viollier, P. H., M. Thanbichle, P. T. McGrath, L. West, 
M. Meewan, H. H. McAdams, and L. Shapiro, 2004, "Rapid 
and sequential movement of individual chromosomal loci to 
specific subcellular locations during bacterial DNA replica- 
tion," Proc Nat Acad Sci 101, 9257-9262. 

Walsby, A. E., 2000, "A square bacterium," Nature 283, 69- 
71. 

Wang, X., X. Liu, C. Possoz, and D. J. Shcrratt, 2006, "The 
two Escherichia coli chromosome arms locate to separate 
cell halves," Genes Dev. 20, 1727-1731. 

Weart, R. B., A. H. Lee, A.-C. Chien, D. P. Haeusser, N. S. 
Hill, and P. A. Levin, 2007, "A Metabolic Sensor Governing 
Cell Size in Bacteria," Cell 130, 335-347. 

Woldringh, C. L., 2002, "The role of co-transcriptional trans- 
lation and protein translocation (transertion) in bacterial 
chromosome segregation," Mol Microbiol 45, 17-29. 

Woldringh, C. L., and T. Odijk, 1999, in Organization of 



1G 



the Prokaryotic Genome, edited by R. L. Charlebois (ASM 
Press, Washington, D.C.), 161-197. 
Wood, W. W., and J. D. Jacobson, 1957, "Preliminary Re- 
sults from a Recalculation of the Monte Carlo Equation of 
State of Hard Spheres," J Chem Phys 27, 12071208. 



